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Abstract 

The decays of n c to K + K~2{n + ir~) and 3(7r + 7r~) are observed for the first time 
using a sample of 5.8 x 10 7 J/ip events collected by the BESII detector. The product 
branching fractions are determined to be B(J/ip — > jn c )-B(r] c — ► K + K~ir + ir~ir + ir~) 
= (1.21 ± 0.32 ± 0.23) x 10~ 4 , B(J/%b -► 7%) • B(t] c -► K*°K*°tt + tt-) = (1.29 ± 
0.43 ± 0.32) x 10~ 4 , and B(J/ip -► 7%) • B(r) c -► Tr+vr-vr+vr-vr+vr-) = (2.59 ± 
0.32 ± 0.48) x 10" 4 . The upper limit for r/ c — > 07r + 7r it + it is also obtained as 
B(J/ip 7%) • B(r] c (pTr+TT-TT+Tr-) < 6.03 x 10~ 5 at the 90% confidence level. 

PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 
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1 Introduction 



The r] c , a 1 S'o state in the charmonium family, was found in the inclusive pho- 
ton spectra from J/ip and i/j(2S) [1] decays, as well as in hadronic decays [2]. 
A number of decay modes of r\ c were then measured [3] . More recent measure- 
ments of hadronic decays of r] c are listed in Ref. [4]. According to Ref. [5], 
the f) c is expected to have numerous decay modes into hadronic final states. 
Although a number of decay modes of r\ c have been measured by different 
experimental collaborations, the number of measured r\ c decay channels are 
few. This means that many decay modes of r\ c are unknown. The 58 million, 
(57.70 ± 2.72) x 10 6 [6], J/ip events taken at BESII provide a chance to ob- 
serve new decays. In this analysis, i] c decaying into K + K~ti + 'k~'k + 'k~ and 
7r + 7T~7r + 7r~7r + 7r~ are studied using J/ip — > 77^. 

The upgraded Beijing Spectrometer detector, located at the Beijing Electron- 
Positron Collider (BEPC), is a large solid-angle magnetic spectrometer which 
is described in detail in Ref. [7]. The momentum of the charged particle is 
determined by a 40-layer cylindrical main drift chamber (MDC) which has a 
momentum resolution of (t p /p=1.78%a/1 + p 2 (p in GeV/c). Particle identi- 
fication is accomplished by specific ionization (dE/dx) measurements in the 
drift chamber and time-of-flight (TOF) information in a barrel-like array of 48 
scintillation counters. The dE/dx resolution is (JdE/dx = 8.0%; the TOF resolu- 
tion for Bhabha events is (Jtof = 180 ps. Radially outside of the time-of-flight 
counters is a 12-radiation-length barrel shower counter (BSC) comprised of 
gas tubes interleaved with lead sheets. The BSC measures the energy and di- 
rection of photons with resolutions of ue/E ~ 21%\^E (E in GeV), = 7.9 
mrad, and a z = 2.3 cm. The iron flux return of the magnet is instrumentd 
with three double layers of counters that are used to identify muons. 

A GEANT3 based Monte Carlo package (SIMBES) with detailed consideration 
of the detector performance is used to obtain the detection efficiency The 
consistency between data and Monte Carlo has been carefully checked in many 
high purity physics channels, and the agreement is reasonable [8]. 



2 Analysis of J /tp — > jrj c , rj c —*■ K + K 77+77 77+77 

These events are observed in the topology 'jK + K~'it + it~it~it~ . Events with six 
good charged tracks and at least one isolated photon are selected. The selection 
criteria for good charged tracks and isolated photons are described in detail in 
Ref. [9] . Each charged track must be well fitted to a helix, originating from the 
interaction region of K xy < 2 cm and \z\ < 20 cm, and have a polar angle 
in the range | cos#| < 0.8. Here R^ is the distance from the beam axis, and z 
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is along the beam axis. Isolated photons are those that have energy deposited 
in the BSC greater than 60 MeV, the angle between the direction at the first 
layer of the BSC and the developing direction of the cluster less than 30°, and 
the angle between photons and any charged tracks larger than 5°. Two of the 
charged tracks should be identified as kaons by combined TOF and dE/dx 
information. 

A four-constraint (4C) kinematic fit is performed under the hypothesis of 
J/ip — > 7ir f i^7r + 7r~7r + 7r~ , and the XyK+K-n+n-n+n- is required to be less 
than 10. To reject background from J/ip — > 'yyK + K~iT + iT~iT + iT~ , X^K+K-Tr+ir-ir+ir- 
is required to be less than xtryK+K-n+ir-w+ir- ■ Background events from J/ip — > 
K+K-^+K-K+K- are eliminated by requiring X* K +K-n+n-n+n- < X 2 k+k-k+*-k+-k- 
and P m iss > 55 MeV/c, where P m i SS is the missing momentum of charged 
tracks. 




o 1 
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m(K + KWjrV) (GeV/c 2 ) 

Fig. 1. The distribution of m K + K - 7T + 7r - 7T + 7T - for selected events. The histogram 
with error bars is from data, the shaded part is the background estimated from 
J/ijj — > anything Monte Carlo simulation, and the curve represents the fitting 
results described in the text. 

After the above selection, the K + K^n + 7Y^n + n~ invariant mass, niK+K-n+w-n+i 
distribution is shown in Fig. 1. A peak at the n c mass is observed. The shaded 
histogram is the background estimated from 58 million J/ib — > anything 
Monte-Carlo events generated with the Lund-charm generator [10]; no promi- 
nent signal in the n c mass region is seen. Also, 100,000 events for the two 
possible background channels J ftp — > K + K~2(n + 7r~) and J/ip — > 73(7r + 7r~) 
are simulated. After final selection, no events remain in the n c mass region. A 
Breit-Wigner folded with a Gaussian to take into account the mass resolution 
of 12.3 MeV/c 2 at the n c and a polynomial background are used in the fit. 
The fit gives 100 ± 26 n c events with a statistical significance of 4.0 a, where 
the mass and width of n c are fixed to the PDG values [11]. 
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Using this sample, we search for the decay mode r] c — > K*°K tt + ti . To 
select K K 7r + 7r~ events, we require that the invariant masses of K + ir~ and 
K-7t + must statisfy \m K% - 0.896 1 < 0.05 GeV/c 2 . After the K*° and K*° 
selection, the K + K~2(-k + 'k~) invariant mass is shown in Fig. 2. A small peak 
at the rj c mass is observed. The background events corresponding to the shaded 
histogram in Fig. 2 are estimated from K*° and K sidebands (0.1 GeV/c 2 
< \m K+w - - 0.896| < 0.15 GeV//c 2 and 0.1 GeV/c 2 < \m K - n+ - 0.896| < 
0.15 GeV//c 2 ), and there is no evident r/ c signal. 45 ± 15 events are obtained 
by fitting the mass spectrum with a Breit-Wigner folded with a Gaussian 
to account for the rj c mass resolution plus a second polynomial background. 
The corresponding mass and width of the r/ c are fixed to PDG values [11]. 
Since the significance of the peak is only 3a, we also give the upper limit for 
r] c — > K*°K* 7i + 7i~. With the Bayes method, the fit of this distribution yields 
65 events at the 90% confidence level. 




2.5 ■^^^■■■li^HH'^ft 

2.75 2.8 2.85 2.9 2.95 3 3.05 



m(K + KW7tV) (GeV/c 2 ) 

Fig. 2. The distribution of m K + K - lv + n - n + n - for rj c — ► K*° K*° it + tt~ candidate 
events. The histogram with error bars is for data, the shaded part is the back- 
ground estimated from K*°(K* ) sidebands, and the curve is the fitting results 
described in the text. 

The J/ip — > 'jK + K~it + it~-7t + it~ sample can also be used to search for r\ c — > 
07r + 7T~7r + 7r~. For selecting a signal, the K + K~ mass, m K + K -, is required to 
be in the region \m K + K - — 1 .02 1 < 0.015 GeV/c 2 . After this selection, no clear 
T) c signal is found in the distribution of m K + K - 7T + 7T - w + 7T - , as shown in Fig. 3. 
Using Bayes method, a fit to r\ c — > K + K^n + n^n + 7Y^ with a Breit-Wigner 
folded with a Gaussian and a polynomial background gives 13.5 r\ c events at 
the 90% confidence level. 

From Monte-carlo simulation, in which the angle (6) between the direction of 
the e + and rj c in the laboratory frame is generated according to a 1 + cos 2 9 dis- 
tribution and uniform phase-space is used for rj c decaying into K + K~2(tt + ty^) 
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Fig. 3. The distribution of m K + K - 7T + n - 7T + n - for rj c 
events. 



/>7r + 7r 7r + 7r candidate 



and (/)2(7r + 7r ), the detection efficiencies of J/ip — > 7^ c (^ c - *■ K + K n + n n + n ), 
J/^> — > 7?7c(?7c - ► K*°K 7r + 7r _ ), and J/^) — > 7^ c (^ c — ► (/>7r + 7r~7r + 7r~) are de- 
termined as (1.43 ± 0.04)%, (1.36 ± 0.04)%, and (1.01 ± 0.02)%, respectively. 
Therefore, the branching fractions obtained are 



;i.21±0.32) x 10" 4 , 



B(J/1> - 


7^c) 


B(J/1> - 




B{J/^ - 


7^c) 


and 




B{J/i> - 


■* IVc) 



K*°K*\ + ir-) = (1.29 ±0.43) x 10~ 4 



K*°K*\+n-) < 1.86 x 10" 4 , 



(f)7l + 7l-7l + 7l-) < 4.72 x 10~ 5 . 



3 Analysis of J/ip — > 7?7 c ,7y c 



7T 7T 7T 7T 7T 7T 



These events are observed in the topology J/^ — > 77r + 7r~7r + 7r~7r + 7r~. Events 
with six good charged tracks and at least one isolated photon are selected. No 
particle identification is required. To suppress background, a 4C kinematic fit 
is performed under the hypothesis 77r + 7r~7r + 7r~7r + 7r~ , and the v 2 is required 
to be less than 10. To reject background from J/ip — > 3(7r + 7r~) and J ftp — > 
3(7T + 7r-)7r°, we require x? 7r+7r - 7r+7r - 7r+7r - to be less than x 2 ^-^-^- and 

Y 2 

A,7r+7T 7T + 7T 7T+7T 7T° ' 

Figure 4 shows the 7r + 7r~7r + 7r~7r + 7r~ invariant mass spectrum after the above 
selection. A clear r/ c peak is observed. The shaded histogram in Fig. 4 cor- 
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2.75 2.8 2.85 2.9 2.95 3 3.05 

m(jrVjrV7iV) (GeV/c 2 ) 

Fig. 4. The distribution of m 7r + 7r - 7r + 7r - 7r + 7r - for selected events. The histogram with 
error bars is data, the shaded part is the background estimated from Monte Carlo 
simulation, and the curve is the fitting result described in the text. 

responds to background estimated from 58 million J /if; — > anything Monte- 
Carlo events generated using the Lund-Charm generator [10], where no i] c 
signal is evident. A fit of the m^+^-^+^-^+T,- distribution, which is shown as 
the solid curve in Fig. 4, using a Breit-Wigner convoluted with a Gaussian to 
represent the signal and a polynomial for the background, yields 427 ± 64 n c 
events with a statistical significance of 6.9a. In the fit, the mass and width of 
n c are again fixed to PDG values [11]. 

The detection efficiency for J /if; — > 7^ c ,^ c — > 7t + 7t~7t + 7t~7t + tt~ is determined 
to be (3.21 ± 0.04)%, by Monte-Carlo simulation with the distribution of 9, 
the angle between the directions of e + and n c in the laboratory frame, being 
generated with a 1 + cos 2 9 and with n c decaying into 3(n + n^) being generated 
with a uniform phase-space distribution. The branching ratio is then found to 
be 

B(J/ip -> 7 ?7 C ) • B(n c -> n+n-n+n-n+n-) = (2.59 ± 0.32) x 10~ 4 . 
4 Systematic errors 

The systematic errors mainly come from the following sources: 

(1) MDC tracking efficiency 
This has been measured with clean channels like J /if; — > AA and ip{2S) — > 
n + n^J/ip, J /if; — > It is found that the Monte Carlo simulation agrees 

with data within 1-2% for each charged track. Therefore, 12% is conservatively 
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taken as the systematic error in the tracking efficiencies for the 6-prong final 
states analyzed here. 

(2) Photon detection efficiency 

This has been studied using different methods with J/ip — > p°7r° events [12]. 
The difference between data and Monte Carlo simulation is less than 2% for 
each photon, and 2% is taken as the systematic error for the photon efficiency 
in this analysis. 

(3) Particle identification (PID) 

This has been studied with J/ip — > K + K~tt°. The efficiency of the PID from 
data is consistent with that from Monte Carlo simulation. The average differ- 
ence is less than 2%. For J/ip — > r yK + K~ n + it~ it + it~ decay, 4% is taken as the 
systematic error from PID. 

(4) Kinematic Fit 

The kinematic fit is useful to reduce background. Using the same method 
for estimating the systematic error as in Ref. [9], the decay mode J/ip — > 
3(7r + 7r~)7T° is also analyzed. The efficiency difference of the kinematic fit for 
data and Monte Carlo is 7.7%. Since the decay of J/tp — > 3(7r + 7r~)7T° is similar 
to the two channels analyzed in this paper, 7.7% is also taken here as the 
systematic error of the kinematic fit. 

(5) T] c parameters 

Although the i] c signal is clear, the number of events is not large enough 
to determine the Breit-Wigner parameters and the background shape well. 
The variation of the fit solution due to changes of the r] c mass and width 
corresponding to the uncertainties in the PDG, as well as changes in the 
fitting mass region used, is taken as a systematic error and listed in Table 1. 

(6) Background 

For i] c — > K + K~7t + 7t~7t + 7t~ , the biggest background comes from 7] c — > K^K^K + K~ . 
When the invariant mass of 7r + 7r~ is required to be within the Kg mass re- 
gion {\m^+ n - — 0.497| < 0.02 GeV/c 2 ), five events remain in the r/ c mass 
region. If all of them are regarded as signal from rj c — > KgKgK + K~ , the 
background from this decay mode is about 5.1%, and this is taken as the sys- 
tematic error associated with background for this channel. No events remain 
for r/c — > K*°K it + it~ and the upper limit is 2.3 events at 90% confidence 
level. Then the uncertainty caused by rj c — > K^K^-k+h^ is 5.1%. 

For the r\ c — > 3(7r + 7r~), Monte Carlo simulation is used to estimate the back- 
ground from i] c — > KgKg7i + n^ . Using the branching fraction for i] c — > K°K 7r + 7r~ , 
obtained from B(i] c — > K + K~tt + tt~) [11], Monte Carlo simulation indicates 
that 33 background events contribute to the r\ c signal. Compared to the 416 
signal events from fitting the mass spectrum, the background fraction is 7.9% 
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which is taken as the background systematic error for this channel. 

(7) Number of J/ip events 
The number of J/ip events is (57.70 ± 2.72) x 10 6 , determined from J/ip in- 
clusive four-prong events. The uncertainty is taken as a systematic error in 
the branching ratio measurement. Table 1 lists the systematic errors from all 
sources, and the total systematic error is the sum of them added in quadrature. 

Table 1 



Systematic error sources and contributions (%) 



Sources 




K*°K*\+tt- 




3(vr+7r-) 


MDC tracking 


12 


12 


12 


12 


Paticle ID 


4 


4 


4 


negligible 


Photon efficiency 


2 


2 


2 


2 


Kinematic fit 


7.7 


7.7 


7.7 


7.7 


7] c parameters 


9.9 


18.6 


14.7 


7.4 


MC statistics 


2.6 


2.9 


2.9 


1.1 


Background uncertainty 


5.1 


5.1 




7.9 


B(<f> -► K+K~) 






1.4 




Number of J ftp events 


4.7 


4.7 


4.7 


4.7 


Total 


19.4 


25.0 


21.7 


18.6 



5 Results 

The decays of r\ c — > K + K^n + n^n + n^ and r\ c — > tx+ n~ tt+ tx~ tt+ n~ are observed 
for the first time, and their decay branching ratios are measured. The upper 
limits of Tj c — > (\>tx+ tx~ tx+ tx~ and r] c — > K*°K* it + it~ are also set at the 90% 
confidence level. To conservatively estimate the upper limit, the systematic 
error is included by lowering the efficiency by one standard deviation. Table 
2 shows the branching ratio results including systematic errors. 

Using the branching fraction of J/x/j — > 7^ c as B(J/ip — > 7?] c ) = (1.3 ± 0.4)% 
from the PDG [11], we obtain: 

B(r] c -> K + K~ 7i + n- 71+71-) = (0.93 ± 0.25 ± 0.34) x 10" 2 
B(r] c -> K*°K*\+7i-) = (0.99 ± 0.33 ± 0.39) x 10" 2 
B(r] c -> K*°K*\+7i-) < 2.36 x 10" 2 
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Tabic 2 



Numbers used in the calculations of branching fractions and upper limits. 



Decay Modes 




e(%) 


Branching Fraction 


J/ip -► JVc Vc K + K-2(Tr + 7r-) 


100±26 


1.43±0.04 


(1.21 ±0.32 ±0.23) x 10~ 4 


J/4> — > 7ry c , ry c — > K"*°.K'* 7r+7r - 


45 ± 15 


1.36 ± 0.04 


(1.29 ±0.43 ±0.32) x 10" 4 


J/V' 7%, % -» K*°ir*°7r+7r- 


< 65 


1.36 ± 0.04 


< 2.46 x 10~ 4 (90% C.L.) 


J/V> -> 7%, r/ c -> 02(7T + 7T-) 


< 13.5 


1.01 ±0.02 


< 6.03 x 10~ 5 (90% C.L.) 


J/V> -> 7%,% 3(7r + vr") 


427±64 


3.21±0.04 


(2.59 ± 0.32 ± 0.48) x 10" 4 



B(r] c -> 07r+7T-7r+7r-) < 5.81 x 1(T 3 (90% C.L.) 



B(r} c -> 7r + 7r-7r + 7r-7r + 7r-) = (1.99 ± 0.25 ± 0.72) x 10" 2 
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